Josephson-vortex-flow terahertz emission in layered high-T c superconducting single 

crystals 
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We report on the successful terahertz emission (0.6~1 THz) that is continuous and tunable in its 
frequency and power, by driving Josephson vortices in resonance with the collective standing Joseph- 
son plasma modes excited in stacked Bi2Sr2CaCu20s+2, intrinsic Josephson junctions. Shapiro-step 
detection was employed to confirm the terahertz-wave emission. Our results provide a strong feasi- 
bility of developing long-sought solid-state terahertz-wave emission devices. 

PACS numbers: 74.72.Hs, 74.50.+r, 74.78.Fk, 85.25.Cp 
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Rapidly increasing applications of the electromagnetic 
(EM) waves in the terahertz frequency range (~0.1-30 
THz) to nondestructive diagnosis, safe medical imag- 
ing, security measures and inspections, high-speed com- 
munications, etc., demand a novel technique of effi- 
cient terahertz-wave generation. Semiconductor EM- 
wave generating elements, covering a wide-range of fre- 
quency spectra, mainly rely on two operation principles: 
sinusoidal current oscillation and the transition of elec- 
tronic states between quantum levels, which are realized 
by electronic and photonic means, respectively. How- 
ever, terahertz-range EM-wave technology has remained 
underdeveloped due to the following reasons: (i) in elec- 
tronics, the transit of charge carriers in a semiconductor 
device inherently takes too long to produce the terahertz 
oscillation, while (ii) in photonics, the terahertz-range 
photon energy is much lower than the room-temperature 
ambient thermal energy. This lack of generating technol- 
ogy in the terahertz band is often dubbed the "terahertz 
gap" . In spite of recent progress in assorted generation 
techniques Q, 0, 0, the terahertz generation remains 
to be a subject of intensive research efforts. 

The present work describes the direct detection of 
the Josephson vortex-flow-induced terahertz-wave emis- 
sion from Bi2Sr2CaCu208+a; (Bi-2212) stacked intrinsic 
Josephson junctions (IJJs). Since the size of the super- 
conducting gap in Bi-2212 is higher than the Joseph- 
son plasma edge both the emission and the detection 
are allowed in a terahertz range In this study, 

the terahertz-wave emission was confirmed by observ- 
ing Shapiro steps in an on-chip stack of IJJs (detector 
stack) as it was irradiated by the emission from the os- 
cillator stack. The frequency of the emitted wave was 
tuned by the bias voltage, while the emission power was 
controlled by the bias current and switching of the res- 
onating modes. 

Naturally grown layered high-T c superconducting Bi- 
2212 single crystals contain IJJs [||, where a 0.3-nm-thick 
superconducting Cu02 bi-layer and a 1.2-nm-thick insu- 
lating BiO-SrO layer alternate along the c axis. Since the 
superconducting electrodes in IJJs are much thinner than 
the c-axis London penetration depth A a fc (^200 nm) the 
phase parameter <f> n in the nth Josephson junction is in- 



ductively coupled to those of neighboring junctions 0,0. 
Thus, the electrodynamics, governed by coupled Sine- 
Gordon differential equations in a system of N stacked 
IJJs, is represented by the same number of standing 
Josephson plasma (SJP) collective eigen modes formed 
along the c-axis direction. Josephson vortices generated 
by a planar external magnetic field tend to move in a tun- 
neling bias current in resonance with the collective SJP 
modes, with the spatial phase distribution ranging from a 
rectangular (in-phase) to a triangular (out-of-phasc) lat- 
tice. The Josephson-vortex-induced emission power can 
be enhanced by tuning the vortex lattice closer to the 
in-phase rectangular mode . 

The conventional solid-state-reaction method was 
adopted to prepare Bi-2212 single crystals, which were 
slightly overdoped. The double-side-cleaving technique 
[lOj was then employed to isolate and sandwich a finite 
number of Bi-2212 IJJs between two Au electrodes, es- 
sentially eliminating the interference from the pedestals 
(stacks of IJJs outside of but coupled to the stack of IJJs 
of interest) and thus allowing the formation of ideal mul- 
tiple SJP modes along the c-axis direction. We fabricated 
two specimens of stacked IJJs (SP1 and SP2) sandwiched 
between two Au electrodes. Detailed fabrication proce- 
dures are described in Ref. H3|. Figs. 1(a) and (b) show 
the schematic of the typical sample configuration and an 
SEM image of SP2, respectively. The oscillator stacks 
were 13.5x1.4 urn 2 [SP1] and 15.2x1.9 /im 2 [SP2], while 
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FIG. 1: (color online) (a) Schematic of sample and measure- 
ment configurations, showing the oscillator stack (left) and 
the detector stack (right) of intrinsic Josephson junctions, (b) 
SEM image of SP2. 
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the detector stacks were 18x1.3 fim 2 [SP1] and 11.6x2.3 
fim 2 [SP2] in their lateral dimensions. 

The lower inset of Fig. 2 shows the oscillation of a 
Josephson vortex-flow resistance (JVFR) of SP1 (iV=35) 
as a function of external magnetic field at T=55 K for a 
bias current of lb— 2 fiA. The field period of the oscilla- 
tion, 0.51 kG, is one half of H [=$ /Ld=1.02 kG]. Here, 
Ho is the magnetic field corresponding to a flux quan- 
tum, <I>o, through each junction of length L (=13.5 fan) 
and thickness d (=1.5 nm). The Hq/2 JVFR oscillation, 
often observed in the S-shaped embedded structures, is 
known to originate from the interaction between a trian- 
gular vortex lattice and the edge potential for the vortex 
entry and exit at the junction boundaries 0,0]. The 
Hq/2 oscillation in our sandwiched oscillator stack thus 
strongly indicates the formation of the triangular lattice 
in its low-bias static limit. 

The upper inset of Fig. 2 exhibits a single vortex-flow 
branch corresponding to the near-static triangular lattice 
in the low-bias region (V<23 mV) and multiple quasipar- 
ticle branches in high biases at T=55 K in a transverse 
field of H—1.2 T; the field denoted by the arrow in the 
lower inset. In this field range, the single vortex-flow 
branch persists down to T—4.2 K (not shown). In a 
higher magnetic field of £7=3.9 T as in the region (ii) 
of the main panel of Fig. 2, however, the single vortex- 
flow branch splits into multiple vortex-flow sub-branches, 
which arise as a result of phase- locking of a Josephson 
vortex lattice on multiple S JP modes 0, 0] . S JP oscil- 
lations in each mode, standing along the c axis, propa- 
gate along junctions with a characteristic collective mode 
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FIG. 2: (color online) /- V characteristics corresponding to the 
collective Josephson vortex states [(i), (ii), and (iii); refer to 
the main text] of SP1 for H =3.9 T at T=4.2 K. The contact 
resistance was numerically subtracted. Upper inset: a sin- 
gle vortex-flow branch (low biases) and quasiparticle branches 
(high biases) in H—1.2 T at T=55 K, corresponding to the 
arrow position in the lower inset. Lower inset: Josephson 
vortex-flow resistance versus external magnetic field of SP1 
at T=55 K in a bias current of 2 fiA. 



velocity 0,0, c„ = cq/\/1 — cos[nn/ (N + 1)]. Here, n 
(=1~N) is the mode index and cq is the Swihart veloc- 
ity, which governs the EM-wave propagation in a single 
Josephson junction. 

In SP1 the JVFR oscillations are limited in the current 
bias range up to h^2 fiA. The multiple vortex- flow sub- 
branches in the bias range of 2 fiA< lb <8 fiA are clearly 
distinguished from the multiple quasiparticle branches 
above i;,~8 fiA, which corresponds to the critical cur- 
rent for the Josephson pair tunneling in the given mag- 
netic field. In general, a stack of IJJs driven by a bias 
current in a planar magnetic field is in one of the three 
different dynamic Josephson-vortex states as illustrated 
in Figs. 2 and 3(a): (i) the near-static triangular vor- 
tex configuration, (ii) the multiple-mode dynamic vortex 
configuration, and (iii) the irregular vortex configuration 
corresponding to the McCumber state dominated by the 
quasiparticle tunneling. 

The maximum vortex-lattice velocity is governed by 
the propagation velocity of a resonating SJP mode in a 
stack of Josephson junctions. Thus, for the bias beyond 
Ib^St fiA, the vortex-flow state transits to the quasiparti- 
cle tunneling state as the flowing vortices pass the state 
of the maximally allowed cut-off voltage, V max , set by 
the highest plasma velocity in a junction. The maxi- 
mum vortex-lattice velocity, c max , iri SP1 was estimated 
with V max to be ~3.8xl0 5 m/s using the relation 0] 
of V max =NHdc max . The value of c max in this sand- 
wiched oscillator stack should represent the fastest prop- 
agation mode among N collective SJP modes. The Swi- 
hart velocity estimated based on the resistively and ca- 
pacitively shunted junction model |l5l Il6| turns out to 
be Co=2.0xl0 4 m/s. Then the corresponding highest in- 
phase mode velocity, ci=3.5xl0 5 m/s, is in remarkable 
agreement with the observed value of c max in SP1. This 
strongly implies that the last sub-branch [the rightmost 
branch in Region (ii)] corresponds to the in-phase rect- 
angular vortex lattice. 

The moving vortex lattice that is in resonance with the 
collective transverse SJP oscillations is theoretically pre- 
dicted to emit highly coherent EM waves at a junction 
edge @, 0. The detection of the Shapiro steps, which 
arise as a result of the resonance between supercurrent 
oscillations and an irradiated microwave in a Josephson 
junction 0], would be the most direct confirmation of 
the electromagnetic emission from the sandwiched oscil- 
lator stack, with information on both the frequency and 
the power of the emission. 

Fig. 3(a) shows the bias conditions of the oscillator 
stack of SP2 (JV=28), Vw=38.9, 51.4, 54.4, and 60.8 mV 
(denoted by filled squares in its vortex-flow branches in 
the dynamic vortex state) and 100.8 mV (denoted in the 
quasiparticle branch in the McCumber state), in a field 
of 4 T at T=4.2 K. In SP2 the detector stack was coupled 
by the bottom Au electrode with the 2-/zm interspacing 
to the oscillator stack. 

The black dldet/dVdet-vs-Vdet curve at the bottom of 
Fig. 3(b) represents the response of the detector stack 
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without any bias to the oscillator stack. The curve cor- 
responds to simple nonlinear I- V characteristics with- 
out sub-branches in the vortex-flow region (not shown). 
The pronounced enhancement of the zero-bias conduc- 
tance is due to Josephson vortex pinning. The vertical 
dotted lines indicate the positions of Vmax of the detec- 
tor stack, where the broad humps are generated by the 
change in the slope of the detector I- V curves at the 
boundary between the vortex-flow and quasiparticle tun- 
neling regions. In order to reduce the vortex entry to 
the detector stack, its width facing the magnetic field 
(2.3 /mi) was designed to be much smaller than that of 
the oscillator stack. In consequence, not all the IJJs 
in the detector stack took part in the vortex flow mo- 
tion. The effective number of vortex-flowing junctions 
in the detector was estimated to be four by comparing 
values of V max between the oscillator and detector stacks 
|19|. Here, we assumed the maximum velocity of the 
vortices was same in both oscillator and detector stacks 
ppf. For a selected fixed bias to the oscillator stack in 
Fig. 3(a), two clear conductance peaks appear in the 
corresponding dldet / dVd e t-vs-Vd e t curve of the detector 
5.3±0.5 (2.8±0), 6.6±0.2 (3.7±0.2), 
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FIG. 3: (color online) (a) Collective Josephson-vortex-flow 
multiple sub-branches and bias voltages V" osc (filled squares) of 
the oscillator stack of SP2 for H=A T at T=4.2 K. The contact 
resistance was numerically subtracted. The two horizontal 
dashed lines divide the I- V characteristics into three regions 
as in Fig. 2. (b) Responses of the detector stack revealed in 
its differential conductance corresponding to the biases of the 
oscillator in (a). Each curve is shifted vertically for clarity. 



7.9±0.3 (4.1±0.2), and 8.5±0.1 (5.2±0.4) mV for four 
bias conditions in order of increasing voltages in the dy- 
namic vortex state, respectively. Here, V^ m and V^t 
are the voltage positions of the primary peak and the sub- 
peak, which are denoted by the upward and the down- 
ward arrows, respectively, in Fig. 3(b). The primary 
conductance peaks are always more evident than the sub- 
peaks. One also notes that the voltage values of peaks in 
each curve shift in proportion to the increase of the bias 
voltage V osc of the oscillator stack. 

To confirm that the observed peaks are indeed the re- 
sponse of the detector stack to the irradiation |2l|, the 
voltages of the conductance peaks (vdet) are plotted in 
Fig. 4(a) as a function of the bias voltage v osc . Here, 
Vosc (= V osc /N osc ) and v de t (= V d et/N de t) are normal- 
ized by the number of junctions involved in the vortex 
flow in the oscillator and the detector stacks (iV osc =28 
and iVd e {=4), respectively. For the primary peaks, the 
values of v osc are in excellent accordance with the val- 
ues of Vdet, which confirms that the primary peaks were 
caused by the Josephson ac response, i.e., integer Shapiro 
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steps I23. Frequencies corresponding to V^ sc =38.9, 51.4, 
54.4, and 60.8 mV are 0.67, 0.87, 0.94, and 1.06 THz, re- 
spectively [2^|. For the bias conditions of V^ sc =54.4 and 
60.8 mV, the humps at V max in Fig. 3(b) smear out, as 
the outward-shifting positions of the conductance peaks 
overlap with those of V max . 

On the other hand, in Fig. 3(b), the voltage positions 
of the conductance sub-peaks turn out to be one half of 
V? et . The relation between the two kinds of peaks in 
terms of reduced parameters is shown in Fig. 4(a). These 
half-integer Shapiro steps are believed to be caused by 
phase locking of the second harmonics of the Josephson 
oscillations on the primary frequencies |24|. For the bias 
voltage of 100.8 mV in the quasiparticle tunneling state of 
the oscillator stack, without the vortex-lattice formation 
along the c axis, no effective vortex-flow-induced emission 
is expected. Thus, as in Fig. 3(b), the response of the 
detector without any Shapiro-step peaks in this bias is 
almost identical to the zero-bias curve represented by the 
black curve. 

The irradiation power, P, is estimated from the ap- 
proximate relation [Hj], P~A7f Vi^l 2 (3 C /2I C , where Aii 
is the height of the Shapiro current step, V\ is the voltage 
value of the primary peak position in iV-stacked junc- 
tions, fl 2 f3 c —(nh/e)(ef 2 /tJ c ), Q=f/f c is the irradiated- 
wave frequency normalized by the characteristic fre- 
quency f c =2eI c R n /h, (3 C is the McCumber parameter, 
e is the dielectric constant, i(=1.2 nm) is the spacing 
between neighboring Cu02 layers, and I c (J c ) is the tun- 
neling critical current (density). The value of Aii is es- 
timated using the relation AJi =AV\ M? \ vr jm, where AVi 
is the voltage width of the primary peak |22| estimated 
at the base level of each curve in Fig. 3(b). Fig. 4(b) 
illustrates that the estimated power irradiated onto the 
detector stack increases rapidly for lower-index modes; 
P=0.69±0.14, 1.01±0.16, 7.5±0.7, and 15.9±0.5 nW for 
the biases of v osc =1.39, 1.84, 1.94, and 2.17 meV, respec- 
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FIG. 4: (color online) (a) Primary (n=l) and sub-harmonic 
(n—1/2) conductance-peak voltages in the detector stack, cor- 
responding to the integer and the half-integer Shapiro steps, 
respectively, with the emission frequencies denoted in the up- 
per horizontal axis, (b) Estimated power irradiated on the 
detector stack for each bias condition in the oscillator stack. 



tively. The emission power from the SJP modes is indeed 
theoretically predicted to grow for lower-index modes 
as the emitted waves become more phase-coherent and 
reaches a maximum for the completely in-phase mode 
corresponding to n=l One also notices that the two 
bias conditions V osc =5lA and 54.4 mV in Fig. 3(a), al- 
though in the same branch, result in a large difference in 
the estimated irradiation power as in Fig. 4(b). This 
takes place because, as the vortex motion approaches 
the resonance condition in a branch, the energy fed by 
the bias is consumed more for an emission rather than 
speeding up the vortex lattice. This observed emission 



from the vortex-flow region well justifies the excitation 
picture of the SJP modes themselves. A more effective 
impedance matching scheme between the oscillator and 
the detector stacks by combining proper antennas in the 
system is expected to generate clearer Shapiro steps with 
higher emission power. 

In the presence of Josephson vortices a Josephson junc- 
tion itself is bound to exhibit a distributed Shapiro- 
step response even to the monochromatic irradiation |'2'2| . 
Thus, the finite Shapiro-step peak width in Fig. 3(b) 
does not necessarily represent the frequency distribution 
in the emitted waves, which can be determined by the 
heterodyne detection instead 0, l^} . 

The vortex-flow terahertz generation in stacked Bi- 
2212 IJJs has the great advantage of continuous and 
frequency-tunable nature of the generated waves. The 
coherency of the terahertz emission from the rectangu- 
lar vortex-lattice configuration is another advantage of 
vortex-flow technique over, for instance, the radiation in 
a quasiparticle-tunneling region where the phase coher- 
ence is lacking between radiations from adjacent junc- 
tions. The direct detection of emission in this study 
provides an additional noble scheme of developing solid- 
state terahertz emission sources |28| that eventually helps 
bridge the terahertz gap. 
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